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We consider Riemannian 4d BF lattice gauge theory, on a triangulation of spacetime. Introducing 
^\ , the simplicity constraints which turn BF theory into simplicial gravity, some geometric quantities of 

Regge calculus, areas, and 3d and 4d dihedral angles, are identified. The parallel transport conditions 

are taken care of to ensure a consistent gluing of simplices. We show that these gluing relations, 

together with the simplicity constraints, contain the constraints of area-angle Regge calculus in a 
5^ . simple way, via the group structure of the underlying BF gauge theory. This provides a precise road 

from constrained BF theory to area-angle Regge calculus. Doing so, a framework combining variables 
^5 ' of lattice BF theory and Regge calculus is built. The action takes a form a la Regge and includes 

the contribution of the Immirzi parameter. In the absence of simplicity constraints, the standard 
7-H 1 spin foam model for BF theory is recovered. Insertions of local observables are investigated, leading 

to Casimir insertions for areas and 6j-symbols for 3d angles. The present formulation is argued to 

be suitable for deriving spin foam models from discrete path integrals. 
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Lattice approaches to describe gravity take several forms. Among them, Regge calculus and spin foams are expected 
| to be intimately related. Both are standardly built with a triangulation of spacetime. The idea of Regge calculus [p 
is to concentrate curvature around [d — 2)-simplices (triangles in four dimensions). In its original from, the geometric 
' variables are the squared lengths of edges. The dihedral angles (the angles between adjacent tetrahedra) are functions 
\& , of the edge lengths, and curvature is measured by the deficit angle around each triangle, which is basically the sum 
C^l ' of the dihedral angles at simplices sharing that triangle. Regge calculus has been revived in the eighties (see Q 
for reviews and references therein) and several alternative formulations have been discussed, in connection with the 
developments of loop quantum gravity and spin foams: first order Regge calculus Q, using triangle areas instead of 
edge lengths (Q for instance), and recently using areas and 3d dihedral angles @. 

As for spin foams, or rather sums over spin foams, they correspond to a way of writing transition amplitudes between 
spin network states in gauge theories @, 0- A spin foam is a colored two-complex joining two spin networks. For a 
single two-complex which may be seen as dual to a triangulation, the sum over spin foams is a state-sum model whose 
data come from the representation theory of some algebraic structure (typically a Lie group). We thus expect to 
derive spin foam models from lattice gauge theories [1(|, exactly as in 3d quantum gravity @] and lattice Yang- Mills 
■ theory 0. This background independent framework is particularly well designed for topological BF theories 

Spin foam models for quantum gravity take advantage of a reformulation of general relativity as a constrained BF 
theory, Plebanski theory [12]. The so-called simplicity constraints ensure that the field B comes from a frame field 

On the lattice, the bivectors, discretizing the field B on triangles, thus contain the geometric quantities, such as 
tetrahedron and 4-simplex volumes and in particular areas and 3d angles. Dittrich and Speziale [f| have argued that 
area-angle Regge calculus is the natural Regge discretization of Plebanski theory. Although appealing, their porposal 
has some practical drawbacks: the expression of the (4d) dihedral angles in terms of the 3d angles takes a complicated 
form, and so do the simplicity constraints. 

We show in this paper that all these expressions are naturally encoded, in a simple way, through the group structure 
of lattice Plebanski theory. This is done by formulating cross-simplicity and the rules for parallelly transporting 
bivectors using only group variables. We then solve the parallel transport conditions for the holonomies in terms 
of the bivectors and additional angles. The latter are recognized as being the dihedral angles and can be seen as a 
discrete u(l) © u(l) connection for an additional gauge symmetry. 

Another interrogation in spin foams and Regge calculus is the role of the Immirzi parameter. More generally, in 
quantum gravity literature, the reader often finds himself in view of different interpretations. While some authors 
bring forward a topological interpretation of the Immirzi parameter (see [f4| . and references therein), its presence 
turns out to be crucial in loop quantum gravity to define the quantum geometry. It also plays an important role 
in the EPRL [l5| and FK7 [l6j], [l7| spin foam models. The recent semi-classical analyses [l8j], [l9| show that it 
disappears from the action at the critical points, in agreement with the continuum. However, as far as we know, it 
has not been introduced in any form of the Regge action. We here fill the first stage of that gap, proposing an action 
(in fact, two, which are quite similar) which includes the Immirzi parameter. Due to the parallel transport relations 
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and cross-simplicity, it is a function of only areas and dihedral angles, and corresponds to the compactified Regge 
action (29| . distorted by the Immirzi parameter. Similarly to the above-mentioned saddle point analyses, the Immirzi 
dependence vanishes on-shell. 

We thus arrive at an improved setting which mixes the advantages of both lattice BF theory and Regge calculus. 
The geometric quantities are clearly identified. The formulation using group elements enables to translate their 
insertions in the path integral as local observables into insertions in the sum over spin foams. For pure SU(2) BF 
theory, we show that insertions of areas and 3d angles respectively translate into insertions of Casimirs on triangles 
and of 6j-symbols on tetrahedra. 

An underlying motivation of this work is to find a framework allowing to derive the known spin foam models from 
the discretized path integral of constrained BF theory. The main approach to build spin foam models for quantum 
gravity relies on the geometric quantization of the tetrahedron classical phase space. This important idea has been 
developped in i2Q] by Barbieri, Baez and Barrett. The EPRL model then imposes the constraints at the quantum 
level, via equations on Casimir operators. The derivation of the FK7 model is different in spirit, but is also based on 



the quantization of the tetrahedron phase space, using coherent states [21] . If these works have given a clear picture of 
cross-simplicity using the geometric quantization of a single tetrahedron, the issue of gluing adjacent tetrahedra is less 
transparent. It would be interesting to directly derive these models from a (constrained) Lagrangian approach. Several 



approaches have been recently developped in that direction. Baratin, Flori and Thiemann 22] have started from the 
Hoist action, thus without need for simplicity, and used a different spacetime discretization. It would certainly be 
interesting to relate their results to the standard framework. Conrady and Freidel [l7| have given a path integral 
representation of the FK7 and EPRL models. It has nice geometric features, which have enabled the semi-classical 
analysis [l^]. It is however, in our view, specifically designed for the new spin foam models, and we would like to 
have a better control of the imposition of the constraints and of the gluing of tetrahedra. Another programme has 
been initiated in [23| , in which a Lagrangian derivation of the Barrett-Crane spin foam model [24\ is given. It shows 
that the main problem of that model is not the implementation of cross-simplicity, but rather the gluing process 
between adjacent tetrahedra. Indeed, the method allowing to recover the BC model is that used to generate spin 
foam models for unconstrained BF-like theories [25| , and its equations of motion fail to yield the expected rules for 
parallelly transporting bivectors in the BC context. The present framework enables to choose and precisely control 
the way these gluing relations are taken into account in the path integral. Since only group elements are used as 
configuration variables (except areas), this formulation is particularly convenient to derive spin foam models using 
Fourier expansions on the groups. 

The paper is organized as follows. In section [JJ we review the standard lattice approach to BF theory, and solve 
the parallel transport relations for the holonomics. In section [Til we extract the geometric quantities of area- angle 
Regge calculus. It is shown that simplicity together with parallel transport naturally contain the expressions of the 
(4d) dihedral angles and of the constraints of area-angle Regge calculus. We also study the action and the role of the 
Immirzi parameter. It is shown in section IIIII that the present setting reproduces the standard spin foam model for 
pure BF theory, and we compute the insertions of areas and 3d angles. 

During the completion of this work, we have learnt that Barrett et al [l9[ have developped a similar analysis, at 
least to extract the deficit angles, but in a quite different context. They are indeed interested in the semi-classical 
behaviour of the EPRL model, which turns out to admit a Regge form. Our point of view is here that the definition 
of the discretized functional integral from constrained BF theory should itself have an interpretation a la Regge, since 
both are basically lattice gravity. It turns out that this idea has also been adopted by Oriti in (26|, from the group 
field theory point of view, with a similar emphasis on parallel transport relations. 

We are all along concerned with Riemannian gravity, so that the relevant structure group is Spin(4). The Lorentzian 
version of the present analysis will be studied elsewhere. 

I. LATTICE BF GAUGE THEORY 

A. The standard discretization 

Let P be a principal G-bundle over a 4d smooth manifold M, spacetime, G — SU(2) or Spin(4). Consider a 
connection A over P, which will be locally seen, as usual, as a 1-form taking values in the Lie algebra g of G (A IJ 
where I, J = 0, 1, 2, 3 are Euclidean 4d indices, and A 1 , i = 1, 2, 3 in the case of SU(2)). The action of the topological 
quantum field theory called BF (27j, [H| is built with a g-valued 2-form field B, transforming under the adjoint 
representation of G: 



S: 



BF 



j Tr(5AF(A)) (1) 
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where F(A) — dA + ^[A, A] is the curvature of A. 

The action is gauge invariant and the equations of motion are: 



d A B = 
F(A) = 



(2) 
(3) 



where d A = d + [A, •] is the covariant derivative. Thus B can be seen as a Lagrange multiplier imposing A to be flat. 
Moreover, the field B can be completely gauged away due to an additional symmetry: B' = B + dA<fi for any q- valued 
1-form <j>, while A is unchanged. The theory has thus no local degrees of freedom. 

To quantize BF theory, we proceed in a standard way by first discretizing the variables. To fill the gap from 
lattice BF theory to Regge calculus via the Plebanski constraints, we choose a setting well-adapted to Regge calculus. 
Consider a simplicial decomposition of spacetime and ask for local frames on 3-simplices, i.e. tetrahedra, and on 4- 
simplices. This means that curvature is concentrated around 2-simplices, i.e. triangles. We will use the dual skeleton 
to the triangulation in an almost systematic way. In the dual picture, triangles are dual to faces, both denoted /, 
tetrahedra to edges, both denoted t and 4-simplices to points denoted v. The boundary of a dual face is made of the 
edges and vertices respectively dual to the tetrahedra and 4-simplices sharing /. The orientations of tetrahedra and 
triangles induce orientations for dual edges and dual faces. 

As far as BF theory is concerned, the case G =Spin(4) is just two copies of the SU(2) case, so that we focus here 
on G = SU(2). The connection is discretized like in usual lattice gauge theory: we consider SU(2) elements g vt which 
allow for parallel transport between local frames. In the dual skeleton picture, the two ends of a dual edge t correspond 
to the two 4-simplices sharing the tetrahedron t. Thus, a dual edge is attached two group elements g v t, one for each 
end. The curvature around a triangle / is thus measured by the oriented product of these group elements all along 
the boundary of the dual face, starting at a base point (reference frame) v*, gf(v*) — Yitcdf 9 vt %'\ ^ v > v> are source 
and target vertices for each dual edge t. The flatness imposed by the e.o.m. then reads: <?/(«) = id 1 for each /. 

The field B, as a 2-form is discretized on the triangles of the triangulation, in a given frame for each of them, bf(v) 
or bf(t) G su(2). Expressions of 6/ in different frames are naturally related by parallel transport along the boundary of 
/. It simply results from the discretisation of the e.o.m. @, integrated along dual edges between different tetrahedra 
or 4-simplices. 

SU(2) gauge transformations change the local frames on tetrahedra and 4-simplices. A gauge transformation k is 
a family of group elements {k t , k v } acting by: 



which is obviously independent of the choice of the base point for each triangle. Because of the simple rules of parallel 
transport for the bf(t)s, one can choose as basic variables the elements g v t and only one bf(t) per face. 

Things are not so simple when introducing the simplicity constraints, and it is then convenient, in order to solve 
them within each tetrahedron, to start with independent bft in the frame of each tetrahedron sharing /. It is as if the 
triangulation were broken up into a disjoint union of tetrahedra. While solving simplicity within tetrahedron, some 
gluing relations are necessary to stick tetrahedra together and carry the information about parallelly transporting b f . 
The idea is thus to use a measure such as: 



k > g vt = k v g vt k t 
k >b f (t) = kt b f (t) k t 



-1 



(4) 
(5) 



The gauge invariant action takes the form of a sum over triangles: 




(6) 



/ 



n n rf3& /* n <%* - ^ sj) 

(*,«) (/,*) (f,v) 



(7) 



in SU(2) matrix notation, and with: 



9tt> = gvt 1 9vt< 



(8) 



1 In fact, due to the use of group elements in the action, and not Lie algebra elements, the action, given below, only catches the projection 
of gj onto the Pauli matrices so that only the sine of the class angle of gj is restricted to be zero. The class angle can thus be or 27r, 
i.e. gf = ±id 
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standing for the parallel transport between two adjacent tetrahedra t and t' through the 4-simplex v, along the 
boundary of the dual face /. dg is the SU(2) Haar measure and the measure d 3 b is to be precised in the following. 
The gluing relations take place at each pair (f,v), often called wedge, since two tetrahedra sharing a triangle / lie in 
a common 4-simplex v, / and v identifying them. 

The structure group being G — SU(2), the norm and the direction of bft have different roles. We will also see that 
the simplicity constraints assign them different roles in order to reconstruct metricity. First, it is clear that the norm 
of bft and that of bft' are equal, Aft = Aft' = Af. it is a SU(2) gauge invariant quantity. Let us use the following 
parametrisation for each bf t as a 2 by 2 matrix: 

b ft = 2 Af ht ■ °, bf t -a = rift cr z nj/ (9) 

with rift € SU(2). a z is the standard Pauli matrix diag(l, —1), and a is the 3-vector whose components are the Pauli 
matrices. Af stands for the norm of each bft, Aj = — 2tr(6j t ), and is clearly independent of any local frame. For 
reasons which will become obvious in the sequel, we will also call it the area of /. Notice that we can equivalcntly 
consider Af e M or R + , since a change of sign of the area can be reabsorbed into the directions bf t without changing 
the action. The direction bf t € S 2 of bft is encoded in the group element rift. However only two parameters of rift 
are relevant since bft is invariant under the right action of the U(l) subgroup generated by a z , rift ~~ > nf t e~^ ea * . 
It is clear in the Euler parametrisation of SU(2), g = e~i aa ' 1 e~^° v e~^ iaz , that the angle 7 is an ambiguity in 
the definition, which should be integrated out. This allows to rewrite the measure using the U(l) and SU(2) Haar 
measures : 

l[dn(A f ) \{d 2 bft Y[S^{bf t -R{gtt')bft') 

f (/,*) (M 

,4* (10) 

Y[dfi(Af) l[ drift J] / Mfv5 S vw(n f } gtt>n f r e^f"') 
/ (/>*) (M Jo 

R simply denotes the vector representation of SU(2). drift is the SU(2) Haar measure and on the right side of this 
equation, it is understood that the Euler angle 7 of each n ft does not play any role since it can be reabsorbed into the 
angles 9f v which are integrated out in the measure. We are not interested in the precise correspondence, including 
Jacobians, between the two above expressions, since we look for the simplest and most natural measure with regard 
to lattice gauge theory, thus using group elements and the corresponding Haar measures as much as possible. The 
remaining ambiguity is the measure dfi(A), which will be fixed to be the Lebesgue measure on 1R by requiring to 
recover the standard spin foam model for SU(2) BF theory (see section [ 



B. BF geometry 

Up to now, the configuration variables are the group elements g v t, the norms Af and the directions bft, i.e. unit 
3- vectors which can be encoded in the group elements rift up to right U(l) multiplication. In particular, the action is 
a function of g v t, Af and bft- The key idea is now to consider the three angles of each rif t as classical configuration 
variables, together with the angles 6f v which are integrated out to define the measure in (fTU|) . In order to extract 
the dihedral angles, we need to solve the parallel transport relations Q for the holonomies gtt> in terms of the 
bivectors. This has been done implicitly in (|I0p . Given an assignment of elements g v t for parallel transport and rift 
for each triangle in each tetrahedron, the gluing condition between adjacent tetrahedra implies that there exist angles 
9f v G [0,47r) such that: 

nft = gwnft>ei e ^ (11) 

which is precisely the content of the delta functions in (fTU|) . As emphasized, the Euler angles 7 of the variables rift 
are unphysical since they do not apear in the expression of bft- However, it seems that they do play a role in (|ll|l . 
The important point is in fact that gw, being equal to: rift e~i efvl7z ril t ,, does not depend on the Euler angles 7. The 
angles 6f v have to take care of the relations between adjacent rift and riff for any local choice of these Euler angles. 

The enlarging of the configuration space thus leads to an additional U(l) gauge symmetry, acting at each pair (/, t), 
on the right olrif t . This means that each pair (/, t) is equipped with a U(l) reference frame which can be arbitrarily 
transformed by U(l) right multiplication. The angles 6f v appear as a sort of connection for this gauge symmetry. 
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Let us note that to take care of the SU(2) gauge invariance, one only has to look at the left of nf t . Indeed, these 
symmetries act on the new variables as: 

(k,X) >n ft = k t n ft ei x ^ (12) 

(M) >0 fv = 9 fv + e{ tl (\ ft -\ ft ,) (13) 

where e£, = ±1. Since t and t' are adjacent, the corresponding dual edges share the vertex v along the boundary of the 
dual face /. e{ t , is positive when the path (t — > t') through v is oriented like /, and else negative. The transformations 
of the angles 9f v are clearly a discretisation of the usual u(l) gauge transformations for a connection. They transform 
so that the holonomies gtt> are indeed independent of the local choices of Euler angles 7. 

Consequently, we will be interested, to make the link with Regge calculus, in SU(2) and U(l) gauge invariant 
quantities. When solving the simplicity constraints, we will show that such angles have the meaning of dihedral 
angles between tetrahedra t and t', and that the U(l) subgroup preserving bft, generated by rift exp(— i0a z /2) njl , 
acquires a real physical meaning. 

To get a hint about the meaning of these angles, note that the content of the delta functions in p0[) can now be 
used to rewrite the action as a function of rift an d 0f v , in addition to and g v t- The relation (jll[) can first be 
inverted to express gw in terms of rift, rift' and 9f v . Thus, it corresponds to (partially) solving the discrete analog 
of the e.o.m. <IaB — @ for the connection, in terms of the discrete B, with a free parameter, 6f v . The holonomy 
around each face is then given by: 

gf{t)=n f te-i<>f°*n- f t =e-^ 6fbft (14) 

where Of = X«e9/ Qfv Thus, Of represents the class angle of the rotation gf(t), and its direction is directly given 

by that of bft, bf t . In particular bft is left invariant by gf(t) as expected. We can now rewrite the action with the 
additional variables: 

S BF (Af,nf t ,0f v ,g vt ) = ^A f sin (^j (15) 

It is trivially SU(2) gauge invariant, and also invariant under the U(l) transformations (I12| . Provided that the angles 
Qf v are interpreted as dihedral angles and Af as the area of /, it corresponds to the so-called compactified Regge 
action, proposed in [29| and (30j . To make this statement precise, we need to study the imposition of the simplicity 
constraints which ensure metricity. 

II. SIMPLICITY AND AREA-ANGLES REGGE CALCULUS 

A. Parametrisation of the constraints 

In the continuum, gravity is described with Spin(4) BF theory when the field B is restricted to a configuration 
space satisfying: 

Bjw Bt L « (16) 

where /z, v, A, a are spacetime indices. When the proportionality coefficient is non- vanishing, the theory consists in 
two branches up to global signs [l3T| : B IJ = ±(e J A e J ) or (-kB) IJ — ±(eJ A e J ) for a non-degenerate cotetrad 1-form 
e. Gravity corresponds to the second case, up to a global sign in front of the action, whereas the first case will be 
called to be non-geometric. Its e.o.m. are just d^e = 0, which constrain the torsion to vanish but leave curvature 
free. 

At the discrete level, the constraints only involve the so-called bivectors Bj J (t), which can be written in terms of 
their self-dual and anti-self-dual components, Bf(t) = b+f(t) © b-f(t) due to the splitting spin(4) = su(2) x su(2). We 
define the Hodge operator * to change the sign of the anti-self-dual sector by the following action on 5pin(4) elements: 
{■kB) IJ = \t IJ kl B KL ■ The constraints are split into three types : 

Diagonal simplicity (*#/(*)) u B 1 / (t) = (17) 

Cross-simplicity (*B f (t)) u B 1 / (t) = (18) 

Closure relation ^e /t 5/(f)=0 (19) 

/cat 
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where eft = ±1 according to the relative orientation of the dual face / and the dual edge t. In the second line, 
/ and /' are triangles of a given tetrahedron. The first constraint, (fT7|) . often called diagonal simplicity, expresses 
the fact that Bf is simple, i.e. it is the anti-symmetrized product of two vectors Bj J — u^v J ^ (and the same for 
its Hodge dual). The second constraint, called cross-simplicity constraint, asks for the sum of two bivectors of a 
tetrahedron to be also simple. Finally, the constraint (jTSJ) imply that all tetrahedra all closed and is thus mentioned 
as the closure relation. When these constraints are satisfied for a tetrahedron, its edges can be labelled by 4-vectors 
only spanning a three-dimensional subspace of 1R 4 . As in the continuum, the bivectors or their Hodge dual are wedge 
(anti-symmetrized) products of these vectors up to a sign~e = ±1. To distinguish between the two solutions, (|18p can 
be changed with: 

eiJKLBy^N? = 0, or B 1 / (i) N tI = (20) 

respectively for the geometric and non-geometric sectors. N t is in both cases interpreted as a unit vector perpendicular 
to the edges of the tetrahedron t. The interesting point is that, in a 4-simplex, simplicity of the five tetrahedra, together 
with the relations of parallel transport between adjacent tetrahedra 

B f (t)=G tt >B f (t')G-}, (21) 

with Gtt> = G^Gvf, induce metricity of the 4-simplex. This is due to the fact that two triangles sharing an edge 
in 4-simplex can always be viewed to be in the boundary of a tetrahedron (this point will be crucial in the present 
work). Equivalently, one can define the bivectors in the frames of the 4-simplices to be: Bf t (v) = eft G v t Bf(t) G~^, 
and use the theorems of [24[ within each 4-simplex and of _18| to reglue them. Notice that the sign e is then a global 
sign. 

Diagonal and cross-simplicity are easily solved by introducing bivectors Bf t — b + f t © &-/t for each tetrahedron 
sharing a given triangle /, and working with their self-dual and anti-self-dual components, parametrized as in the 
previous section ([9]), b±f t = l>A±f n±ftcr z n^ t . The areas A±f do not depend on tetrahedron frames since the Spin(4) 
parallel transport does not affect the norms. Then, ([!?)) and (|18l) can be solved independently within each tetrahedron, 
implying relations between the self-dual and anti-self-dual areas A±f and directions n±f t . Indeed, (|17p means that 
b±ft have equal squared norms, A 2 _^ = A 2 + j. Equivalently, they are related by a rotation hft £ SU(2): 

3/i/ t eSU(2) b-ft = ±h f - t 1 b +f thft (22) 

As shown in [2l[, in the geometric (respectively non-geometric) sector determined by the sign '-' (resp. '+'), hf t has 
a clear geometric meaning. The Spin(4) element Hf t = (hft, id) maps the reference vector A^°) = (1,0,0,0) to a unit 
vector Nf t orthogonal to / in the sense (*B ft ) IJ 'Nfti = in the frame of t (resp. B^Nfti = 0). 

hft thus stands for the choice of a normal to /. Cross-simplicity then imposes that such a choice can be done at 
the level of each tetrahedron, for its four triangles alltogether fl6j.pH: 

b„ ft = ±h- l b + fth u (23) 

h t representing a unit vector AT/ = (h t , id) 1 J Nj > ^ satisfying the constraint (|2TJ)) . A gauge transformation K acts 
independently on the self-dual and anti-self-dual sectors according to It acts on h t so as to preserve the Spin(4) 
covariance while imposing the constraints (which relate the self-dual and anti-self-dual sectors) : 

K > h t = k+ t h t kz] (24) 

Due to this specific transformation, h t can always be gauged-fixed to the identity, which corresponds to the usual 
time gauge in loop quantum gravity. In this gauge, all tetrahedra are orthogonal to the reference vector . This 
does not imply trivial correlations, since quantities must be compared in a common frame. For instance, comparing 
the normals of two adjacent tetrahedra t and t' can be done by transporting h t > in the frame of t, i.e. g+tt' fh' 9-tt n 
After the gauge-fixing, correlations are encoded into g±tr ■ 

Relations (|22p and (|2"3"|) imply that hft and h t differ by an element of the U(l) subgroup leaving b + f t invariant. 
Indeed, a triangle in 4d admits an orthogonal plane on which hf t and h t represent different choices of unit vectors. 
This orthogonal plane is generated from a normal of reference by action of the U(l) subgroup preserving 6+/t- 



2 To eliminate this sign ambiguity, one should ask for a condition on the positivity of tetrahedron volumes, as in |24l| . We will not deal 
with such a constraint, as it is not present in the known spin foam models for quantum gravity. 
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Consequently, as a triangle is shared by several tetrahedra, the vectors N t are not independent of each other. For 
two adjacent tetrahedra t and if, N t and N t < both lie on the plane orthogonal to *Bf t (in the geometric sector), up 
to parallel transport, provided regluing conditions hold so that the bivectors Bf t and Bf t i are related following (12ip . 
This geometric information is recorded through an angle ipf v G [0,47r) and the relation: 

fh (g+tt< h t > gZ],) ~* = n ft nj? (25) 

Having introduced the normals Nt, the constraints can be formulated so as to be linear in the bivectors (|20 |) . (|23|) . 
Without an explicit regluing such that presented in the previous section, equation (|23[) sets up the framework of the 
Barrett-Crane spin foam model. The linearity of the action and the constraints in the bivectors together with the 
formulation of the normals Nt as group elements make the computation of the partition function easy to perform. The 
regluing is implicitly (and incorrectly) performed via boundary variables for the discretised connection (see details in 

Before gathering the gluing relations of the previous section with the simplicity constraints, we need to rewrite the 
latter in terms of the area and direction variables. We have self-dual and anti-self-dual areas satisfying A + j = ±A_/, 
self-dual and anti-self-dual rotations n±f t encoding the directions of 6±/t, which are related by the group elements 
h t . Taking into account the U(l) x U(l) invariance in the definition of Bf t — (n+/t,n_/t) : 

Hd»(A +f )dn(A^ f )6(A +f TA- f ) 11 dn +ft dn_ ft f J[dh t J] / #/* S ( n -)t K 1 n+ ft e****") (26) 

/ (/,*) * (/,tr° 

Diagonal simplicity and the choice of the sector are contained in the one-dimensional delta functions on the areas. 
Cross-simplicity is imposed by the last SU(2) delta functions in (f2"6"|) , and dh t is the SU(2) Haar measure. The closure 
relation has to be added. Due to cross-simplicity, only its self-dual components are needed : 

II S { £ /* A +f n +ft^ n +}t) ( 27 ) 
t fedt 

Taking two copies of (|10[) the regluing relations imposing ([21]) read : 

/>47T />4?r 

n / de+j(nl) t g +w n +ft ,ei e ^) / dSj v 8{nZ)t9-w n. fv e^*) (28) 

(f,v) J ° J ° 

They obviously imply that Bf t is left invariant by the holonomy based at t, Gf(t) Bj t Gj (t) = Bf t . As a consequence, 
the normal vector Nt is not preserved by Gf(t) but send to another which is still orthogonal to *Bf t : 

(*B ft ) I J (G f (t) J K N t K )=0 (29) 



B. Geometry and relation to area-angle Regge calculus 

After the simplicity constraints are taken into account, the initial variables Bf t and G v t = (g+vti 9—vt) contain all 
the geometric information. It is then interesting to know how the transition to a certain form of Regge calculus can be 
explicitly made. We may suspect the angles 9 ^ v and ipft introduced in the measure to be related to dihedral angles. 
Moreover, having partially solved for the discretised connection with the gluing relations, we expect to switch from 
the first order setting of BF theory to a second order formalism, provided 4-simplices are not degenerate. 

The variables of area-angle Regge calculus are the triangle areas and the 3d dihedral angles 4*^, , i.e. the angles 
between two triangles /, /' within the tetrahedron t. Areas and 3d angles are naturally subject to constraints, which 
are going to be discussed below. The 4d dihedral angles U ' , measuring the angles between tetrahedra and which 
appear in the Regge action, are in [5J functions of the 3d angles. Using the geometric meaning of the variables Bf t , 
h t and G vt , we are going to define some quantities which will naturally be interpreted as areas, 3d and 4d dihedral 
angles. We will then check that the gluing relations together with simplicity lead to the standard expected relations 
between these quantities in area-angle Regge calculus. We emphasize that the relations derived in this section are not 
new complicated relations to be taken care of the path integral, but are simply consequences of the parallel transport 
conditions (|2"8|) and of the simplicity constraints as written in (|2"6"]) . 

Like in the previous section, instead of the true directions &±/t of the bivectors, which are in S 2 x S 2 , we consider 
the configuration variables as being n-ft) G Spin(4), together with the angles 9^ v and ipft which are integrated 
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in ([2"o]) and (f2"5)) . This enlarging of the configuration space enables to solve the parallel transport relations for the 
holonomies and to express cross-simplicity only in terms of group elements. That was already the content of the 
delta functions in the above measures. Given an assignment of G v t and (n_|_ n_/t), the regluing and the simplicity 
constraints ask for the existence of angles such that : 

g+tf = n +ft e _ 5 e /^ n -) t , : g-tf = n_/ t e"* /-^ nZ) t , and n_ /t = n+/ t e5*"^ (30) 

Let us emphasize that n±ft and the angles do not catch all the information about the discrete holonomies, since Gtt' 
is defined as Gtt' = G~l G v y ■ Moreover, the group element G v t is the same for the four dual faces sharing the dual 
edge t, while n±ft and the angles depend on those faces. 

Like in pure BF theory, the Euler angles 7 of n±f t are unphysical (for n — e% aaz e% ,3 ' Jy e^ 7 "^) . The crucial point is 
that g± U ' and h t , as expressed in (|30p . must not depend on local choices of those Euler angles. The enlarging of the 
configuration space is thus compensated by an additional U(l) x U(l) gauge symmetry, very similar to (H2): 

(if,A) >n ±/t = fc ±t n ±/t e* A /* CT * (31) 
(K, A) > 9% = 0% + 4, (X% - X%,) (32) 
(K, A) > ipf t = ijjft - (X]~ t ~ Xj t ) (33) 

The angles 0j v and 9J V form a sort of discrete u(l) © u(l) connection. The angle ipft transforms so as to preserve 
the U(l) x U(l) covariance while solving cross-simplicity. Its role is similar to that played by the rotation ht for the 
SU(2) x SU(2) gauge symmetry, iffij]). In particular, we may choose a gauge fixing such that ht = id and tpf t = 0, 
so that the cross-simplicity constraint, i.e. the third equation of (|30|) . reduces to 8{n_ft n +ft)- We now consider the 
geometric quantities of interest, which are both Spin(4) and U(l) x U(l) invariant. 

The simplest local observable is the area of a triangle. Since in the continuum the field B is built from the cotetrad, 
in the lattice the norm of squared Bf t represents the squared element of the metric integrated in the directions of the 
triangle, i.e. its area. Because of diagonal simplicity, A 2 +f = A?_f, the self-dual part is sufficient to recover the area: 

TT(B%)=2tr(b +ft )=-A\ f (34) 

Considering the products of bivectors within a single tetrahedron then enables to reconstruct its geometry. Indeed, 
the angle <p 1 ^, between / and /' in t is related to the dot product of the 4-vectors N t j B 1 ^ which are orthogonal 
to triangles in the 3d subspace spanned by t. It is easy to see that those dot products correspond to looking at the 
traces: TrBf t Bfi t [Hj]. Because of cross-simplicity, the self-dual components are sufficient. We precisely consider: 

cos <p) f , = -e ff , b +f t ■ b+f> t (35) 

= D oo ( n + f t n +f't) ( 36 ) 

where £//' = £/*£/'* is the relative orientation of the dual faces / and /'. The last equality is a simple computation 
involving two spins 1/2 which couple only in the spin 1 representation, b ■ b' = ^tr(na z Ti~ 1 n'er 2 n' -1 ) = D^ (n _1 n'), 
where Dq 1 ^ is the diagonal matrix element in the spin 1 representation for the state |1,0). Using the Euler angles, 
g = e~i aa * e~iP av e _ ^ 7<T * , it is a very simple function which is in fact the Legendre polynomials P\. (l,0|.g|l,0) = 
Pi (cos 0) = cos p. As expected, the U(l) action on the right of each n + f t is absorbed by the state |1,0). The 3d 
angle (pff l& thus a SU(2) and U(l) gauge invariant quantity which is encoded in the correlations between the group 
elements n+ft and n+fn- 

As for the internal (4d) dihedral angles between two adjacent tetrahedra, they are defined by the product of the 
normal vectors, — Nt ■ Nt> , up to parallel transport: 

cos0 tt , =-N t (v)-N t ,(v) (37) 

= -i tr h t g-tt> h", 1 g~} v (38) 

The last equality simply follows from the standard action of Spin(4) on R 4 . The key quantity to be studied is thus 
the group element H t t', defined in the frame of t, which compares the two normals: 

H w = h t g-w K, 1 gl\ v (39) 



9 




f 



FIG. 1: In the dual picture, the three tetrahedra t, t' and t" of the 4-simplex v become edges meeting at the vertex v. Each 
triangle being shared by two tetrahedra in v, the boundaries of the dual faces in the neibourhood of v are made of two dual 
edges. For a pair of tetrahedra meeting at /, one can equivalently write the dihedral angle 8 tt n between them in terms of the 
3d angles using three different intermediate tetrahedra (the boundary of a 4-simplex being made of five tetrahedra). This leads 
to the constraints (|48[) between the 3d angles, as proposed in [f| 



As we have already seen, it can be expressed as an element of the U(l) subgroup which leaves 6+/t invariant (1251) . 
since N t and N t i are both in the plane orthogonal to *Bf t . We are thus mainly interested in the class angle. Using 
([50)1 , it can be expressed in terms of the angles 9^ v and ip ft ■ 



+ft 



(40) 



where t{ t , = ±1 is positive when the path (t — » t') through v is oriented like /, and else negative. Let us define the 
internal dihedral angle 9 t t>, whose label is equivalent to 0f v , via: 



cos 6tt> = — cos 



sm Vtt> 



- C J 



e ey sm ; 



(41) 
(42) 



where e is precisely the global sign ambiguity appearing in the expression of the bivectors in terms of the edge vectors 
[l||. It is then convenient to express Hu> with the direction and the dihedral angle between t and t': 



cos 



it sm 



(9tf)'b 



+ft 



(43) 



with b+ft ■ 5 = n + f t <7 z n +ft . 

We have started with lattice BF gauge theory, which is a first order formalism, in which the connection is independent 
of the field B. In the continuum, when simplicity is implemented, the equations of motion cIaB = reduces to 
d,A& = if e is a non-degenerate cotetrad. This latter equation then admits a unique solution A(e), which is said 
to be compatible with the cotetrad. Here, the gluing conditions have partially solved a discrete analog of the e.o.m. 



cIaB = for the discrete connection, g±t 



i±fte 



L ±ft' 



Barrett [3(] has emphasized that first order Regge 



calculus can be denned provided the angles Qtv are restricted to be dihedral angles of any well-defined geometry of 
a 4-simplex. Since here the only geometry is that given by the bivectors, the natural question for our setting is to 
know if the simplicity constraints together with the gluing relations are sufficient to relate the dihedral angles 9w to 
the geometry descibed by the bivectors Bft, i.e. to the 3d angles ^ffi according to (|44"|) . This is indeed the case. 

Consider three tetrahedra t, t' and t" in a 4-simplex v, each pair meeting at a triangle, according to fig. [T] The 4d 
dihedral angles are defined in a flat 4-simplex in terms of the 3d angles by: 



cos 0* 



sm 



sm 



(44) 



It is important to see that the intermediate tetrahedron t' , with triangles f\ and /2, can be changed with another 
tetrahedron of the 4-simplex. This leads to constraints for area-angle Regge calculus between the 3d angles. The 
complicated relation ()44|1 is in fact naturally encoded in the present framework into the constraints written as relations 
among SU(2) variables. To obtain it, one can recognize the structure of the multiplication law of SU(2). The idea is 
then to express H U " in terms of H U ' and H t >t"- From its definition (|39p . one simply has: 

H tt n = H tt , g +tt , H VV > g+] t , (45) 

This is the key identity in our setting. From the definition of H t f, it is quite trivial and mainly due to the fact 
that parallel tansport between t and t" through their common triangle / can be performed using t' as a pivot, along 
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the boundary of the dual faces /i and / 2 instead: g +tt " = g + l t g +v f g + l t > g v t" = g+tr g+t't"- However, when H tt ' is 
expressed as a U(l) element preserving (|43p . it becomes a non-trivial statement. Taking the trace of equation 

(US]), using the regluing conditions along /i between t and t' , one obtains: 



cosfl tt / = -tr 



' ' cos6» tt - +iee(l, sm6 tt > n + / lt /<7 x n + ^ t/ ) (- cos^/ t » + ef, 2 t „ sin6» t - t » "+/ 2 t'^^ + } 2t /) (46) 



from which it follows that: 

cos tt 



>=-cos6 tt ' cos 6 t 't" + sin 9 U ' sm9 t ' t » cos<^ l/a (47) 

This relation is naturally SU(2) and U(l) invariant. For non-degenerate configurations such that sm9 tt > ^ 0, it gives 
the 3d angles as functions of the 4d angles. It can also be inverted, precisely yielding to (|44|) . Interestingly, since (|45|) 
is an equation between SU(2) elements, its projection onto a also gives relations between the 3d and 4d angles, which 
are only SU(2) covariant. However, their geometric meaning is unclear to us. 

We have almost completed the proof that (|30"|) and the closure relation contain the constraints of area-angle Regge 
calculus. Indeed, the expression of the angles 8 tl t 21 (]44[) . has been established for any intermediate tetrahedron. The 
3d angles hence satisfy: 



cos ^ h - cos 4>)\ f cos ct>% _ cos 4 f , - cos / cos 4*, 



sin 4>)\ f sin 4 2 /2 sin 4>)\ f sin 



(48) 



which implies a consistent gluing of tetrahedra. Here, t\ and ti share the triangle /, t\ and t share /i, and t and 
ti share ]i (and similarly for the primed labels). Moreover, by contracting the closure relation (|19[) with a specific 
bivector ef t Bf t , we obtain a gauge invariant form: 

A+ f - A +f 008 rtr =° ( 49 ) 

f'^fedt 

The constraints (|4"B")) and ([4^]) are precisely proposed in [5| as a set of constraints for area-angle calculus. 



C. The action and the Immirzi parameter 



Before studying how these variables can be useful to build some old and new spin foam models, we need to write 
an action. It is expected to look like a (compactified) Regge action, but the role of the Immirzi parameter in this 
context is badly known (in particular, it was not introduced in Q). Spin(4) BF theory can be equivalently formulated 
with the field B or its Hodge dual. It turns out that even for gravity, when B — *(e A e), a term proportional to 
tr(*.B A F) can be added to the action without modifying the e.o.m. (at least for pure gravity, without matter), the 
proportionality coefficient being known as the inverse of the Immirzi parameter 7. Defining j± = ±(1±7 _1 ), the 
naive action is thus: 

/- i - = ^tr((l + 7 - 1 *)% i G / (i)) =£ 7+ tr(6 +/tff+/ (i)) - 7 -tr(6_ /t £_/(<)) (50) 
/ / 

For each triangle, this action needs a base point, i.e. a tetrahedron of reference. The gluing relations between bivectors 
of the same triangle, taking into account the structure of parallel transport, makes it independent from the choice of 
the base points. However, we clearly see that simplicity, especially when rewritten as in (|2"3")l . leads to traces of 6+/t 
times a linear combination of SU(2) elements. This unnatural situation can lead to spurious measure factors in the 
functional integral, as shown in 23]. It is far more natural to use instead the group structure, which corresponds at 
first order to the sum of the Lie algebra elements. In this context, the Immirzi parameter is introduced as a coefficient 
for these algebra elements, as already proposed in [23j. Let us define, for g = cos 9/2 + i sin#/2 h ■ a € SU(2) of class 
angle 9 and direction n, the element of class angle a.9 and same direction: 

9 9 

g a = cos(a-) 4- i sin(a-) n ■ S (51) 

For integral a, it is obviously the group multiplication. To introduce the Immirzi parameter, we may propose: 
2 f tr(6+/t g+t' 1 it)) + tr(&_/t g_ * 1 (t)). The cross-simplicity constraints (|23|). b-ft = —h^b+ftht, still lead 
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to sums of group elements: tr b+ft (t) — htg_J (t)h t ). We thus change that expression by using the group 
multiplication which corresponds at first order to the sum of Lie algebra elements like in the continuum action. 

I 7 (b +fu ht,g± vt ) =J2 t *( b +ft9Vf(t)h t gl- f (t)ht 1 ) (52) 
/ 

The regluing relations and the simplicity constraints ensure that it is gauge invariant and independent from the choice 
of the tetrahedron of reference for each triangle. Note that ([52")) uses cross-simplicity for each tetrahedron of reference. 

The sign '-' in front of (1 — 7 _1 ) for the anti-self-dual g_ ^ 7 ' is due to the sign '-' in the constraints for the geometric 

case, 6-/t = —h^ 1 b+ft ht- Its removing corresponds to the choice b-ft = h^ 1 b+ft ht- 

Another proposal is to impose cross-simplicity for each tetrahedron of reference via the group multiplication, like 
in (|52[) . but to keep separately the contribution of the Immirzi parameter: 

Iy(b + ft,h u g±vt) = ^tr(b +f t g+ f (t) h t gZ){t) K 1 ) + 7~ 1 tr(& +/f g +f {t) h t g-f(t) h^ 1 ) (53) 
/ 

It is also SU(2) gauge invariant and independent from the choice of the base points. 

Notice that the naive continuum limits of 7 7 and J 7 are exactly the same as that of the naive action J" 3,1 ™, (|50|) . 
Consider a chart where the typical length of edges of / is of order e. When e goes to zero, we can use the expansion: 
g±f w 1 + e 2 F±\f, where F±\f is the component of the curvature along the directions of the face dual to the triangle. 

This reduces /" alvc and I 7 to the same expression: J2f 1+ £2 ^ jT Q ) +ftF+\f(t)) — 7-£ 2 tr(&-/ti* 1 -|/(£)). Also, with the 
Immirzi parameter, g 1 ^ w 1 + 7±e 2 -F±|/- Inserting this expansion into (|52p . we are lead to the same result to first 
order in e 2 , and thus to the same continuum limit. 

If we again promote the angles 6f v and tpft appearing in the measure (128|) . (|26[) to configuration variables, we 
can easily express the action 7 7 with only areas and angles. Indeed, from (|30[) . the holonomies around triangles are 
determined by the directions of the bivectors and the self-dual and anti-self-dual deficit angles 9^ = J2v&af ®% : 

g +f (t) = n +ft e -i 9 l a > n-) t and g_ f (t) = n_ /t e^ 9 !^ n~_) t (54) 
The angles of the rotations g 1 ^ are simply 7±0y • This leads to the following simple forms: 




and for / 7 : 

J 1 {A + f,k + f U e%,ht,^ fu g± v t) =J2 A +f sin ^/ - J) +^ A +f sm^(9+ + 9j) (57) 

/ 

As expected, A + f represents the area of the triangle / as a function of the edge lengths once simplicity is taken into 
account. We recognize the usual, geometric, deficit angles, Of = (6j — 9j)/2. Notice that Of is in fine insensitive 
to the angles ipf t appearing in the simplicity constraints (|3H)) and which play a role in the dihedral angles (|4*01) . But 
we also see the appearance of some non-geometric deficit angles, given by the sum of the self-dual and anti-self-dual 
angles, instead of their difference, due to the presence of the Immirzi parameter. 

In the continuum, the e.o.m. <1aB = with B = *(e A e) for a non-degenerate cotetrad leads to rf^e = 0. In 
the lattice and for nondegenerate sets of bivectors, simplicity and parallel transport make sure that the edge vectors 
ei(t) defined in local frames of different tetrahedra are related by parallel transport, up to a change of sign [l8[. In 
particular, the holonomy G/ around a dual plaquette leaves the edges of the corresponding triangle invariant up to 
signs: 

Gf{tYjez{t) J = e i ^^e l {t) 1 (58) 

where r\t = 0, 1 and the sum is over the dual edges around the dual face /. Again in the continuum and for a 
nondegenerate cotetrad, the equation d^e = can be solved giving a unique connection A(e) compatible with e. This 
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connection is such that the Immirzi term of the action, e 1 A e J A F]j(A), identically vanishes. We would like our 
discrete setting to also consistently imply the disappearance of the Immirzi parameter from the e.o.m. (involving 
constraints between the angles 0~j and 6j) for nondegenerate configurations. The parameter 7 only appears in the 
action for the closed holonomies Gf(t), so we do not need to carefully look at each step of parallel transport between 
adjacent tetrahedra. Without taking (|58p into account, Gf(t) only leaves Bf t invariant and contains a priori terms 
proportional to (0j + 8 J ) as well as Of — (9~t — 0j)/2: 

G f {t) = e^P+f'- 6 ^-?^ (59) 
= e 2^ e /+ 9 7)* s /' e-xti 6 *** (60) 



When simplicity is imposed, the edge geometry and in particular equation (|58| have to be considered. It is clear that 
the terms proportional Of are not involved, since Bj^ egj(t) — in the geometric sector. Thus, the non-geometric 
angle (Oj + 9j) is responsible for the phase factor of (|58|) : 

0+ + Oj = 2vr Vt (mod 4tt) (61) 
t 

Gfit) = e A +f * * ft e A +f 1 }t for non-degenerate configurations. (62) 

When the phase factor is trivial, this is a well-known property of Regge calculus. Considering a plane in R 4 generated 
by the vectors e\ and e2, only a U(l) subgroup leaves any vector of the plane invariant, which is that generated by 
B\2 = *(ei A e%). The plane is stable under the action of the rotations generated by *B\i = e\ A 62 and (f6"2")l allows a 
reflection, ei.2 — > —61,2, i-e. a rotation of angle 7r generated by *-Bi2- The form of the on-shell holonomy (|62[) implies: 

K 1 g +f (t)h t g- f (t) = e i ^nd (63) 

This can be interpreted in the time gauge, h t = id, in which the edge vectors of t span the 3d space orthogonal to 
For a triangle /, the holonomy takes the form Gf(t) = (g+f(t),g^f(t)), or G/(t) = (g+f(t), — g^\(t)). In the 
first case, it is a boost, in analogy with the Lorentzian terminology, whose axis is the normal Nj to / within the 3d 
space of t, i.e. a rotation generated by 7V( ) A Nj. In the second case, it is combined with a reflection of /. 

Equation (|63[) has also interesting consequences in the dual picture in terms of normals to tetrahedra. Given the 
rules for parallelly transporting h tl (|24p . this equation means that the normal vector N t is left invariant by the 
group element -kGf(t) = (g+f(t),±gZf~{t)). A similar relation was obtained in [23], while dealing with Spin(4) BF 
theory constrained by diagonal simplicity only, in the non-geometric sector. In this situation, h t is replaced with hf 
which stands for the choice of a normal Nf = (hf ,id)N^ to the triangle in the sense Bj J Nfj = 0. The equation 
hj (t) g+f(t) hf(t) g~f(t) = id then results from varying the action with respect to the bivector. In the geometric 
sector, one has instead hj 1 ^) g + f{t) hf(t) g~Ai) — id, meaning that Nf is left invariant by G f. It is thus natural that 
(|63[) entails the disappearance of the Immirzi parameter of the action, when seen as an e.o.m. of the non-geometric 
sector. Notice indeed that the quantity of the l.h.s. of (f63|) is precisely that entering the action 7 7 , (|53|) . with the 
coefficient 7 _1 . It thus becomes on-shell independent of 7: 

I 7 = A +f sinO f (64) 
/ 

As for the action J 7 , it looses its dependence on 7 and reduces to (|64[) when 7 _1 is an integer. 



III. SPIN FOAMS FOR BF THEORY AND INSERTIONS OF LOCAL OBSERVABLES 



In the previous sections, we have discussed an improved setting to discretize gravity from BF theory with constraints. 
We have proposed to explicitly introduce the parallel transport relations for bivectors in the discretized functional 
integral. We here check that in the absence of constraints, the standard spin foam model for BF theory is recovered, 
fixing in addition the integration measure for the area. Moreover, the present framework allows to insert local bivector 
observables, corresponding to areas and 3d dihedral angles when simplicity is fulfilled, and explicitly compute their 
spin foam quantisation within SU(2) BF theory. It is certainly not the most straightforward road to obtain those 
results, but it introduces the methods to be used in the presence of the simplicity constraints and which then turn 
out to allow a straightforward derivation the new spin foam models [38[ . 
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A. BF theory 

As we have seen, the field B in SU(2) BF theory can be seen as a multiplier imposing the vanishing of the curvature. 
Thus, the integration over B in the partition function projects onto these configurations. At the discrete level, the 
holonomy around a dual face / is constrained to be ±id, or equivalently, its projection onto SO (3) is the identity: 



^BF 



/ II dg vt Y[S so(3) (g f ) (65) 
(t,v) f 

= E / U d ^ EK*(S/) ( 66 ) 

{j f GK} (t,o) / 

where dg is the SU(2) Haar measure, and dj — (2j + 1) and Xjid) are respectively the dimension of the representation 
of spin j and its character. If g — e l ^ n ' a for 4> £ [0, 27r), then: Xjis) — b ' s "n • 

Using the variables introduced in the previous sections, the partition function is defined by: 

Zbf = f ]Jdg vt l[dn ft JJ ^(n^ 1 n /t * e i e /^)] JJ d/x^e^ 81 " 1 ^ (67) 
(*,«) (/,«) C/.«) / 

where dn is the SU(2) Haar measure, dO the Haar measure over U(l) (between and Att), while d^,{A) will be chosen 
so that this expression reproduces the standard result ([65]). 

It is clear that for each face the delta functions can be used to eliminate all rift but one, whose tetrahedron is the 
base point, say U/j*, giving: S(nft* <?/(£*) n ~ft* ^ e}IJz ). Integrations over the angles 0f v can be reduced to an integral 
over a single angle Of. In analogy with Regge calculus, it can be said to represent the deficit angle since it is the class 
angle of the holonomy gf(t). 

The spin foam formalism is a way to write down the partition function as a state-sum whose data come from the 
representation theory of the considered groups. We thus need to expand the ingredients of Z BF into U(l) and SU(2) 
harmonic modes. First, for each face, the regluing conditions lead to a delta over the conjugacy class 6/ for the 
holonomy: 



dn 5{ngn- l e^) = T Xj (g) Xj (e^"*) (68) 

SU(2) M 

J 2 

In particular, Xj(^ eaz ) = 5p=-j e ~ w6 ■ 

As mentioned in section |Tf the area can be taken in 1R or in R + . If Aj is integrated on R with the Lebesgue 
measure, then Of is constrained to be or 2ir. It directly follows that gj = ±id and so (|65|) . It is however instructive, 
and in the spin foam spirit, to use expansions whose data can be interpreted and considered as physically relevant. 
Keeping in mind the issue of including the simplicity constraints, we consider Af £ R + and expand the exponential 
of i times the action over U(l) modes: 

AAsin % \ ^ t / a\ Am 



e 



m£TL 



Jm(A)e m i (69) 



This is the Jacobi- Anger expansion, which defines the Bessel functions J m of the first kind: J m (A) = 
Jq^L: cos (mO ) e lA cose . The Bessel functions satisfy some interesting properties [3l|. Their integral over R + is 
normalized to 1 for m £ N, and important symmetries are: J_ TO (A) = J m (—A) = (—l) m J m (A). Combining those 
results, we have: 

dA el ^si„ | = 1 + ( e '»"f + e - im %) (70) 



R + m>0 

Integrating then leads to, for j £ N/2: 

f4w 



EE/^ e-^(e""i + (-ire-^l) = EE + 

m>0 p=—j m>0 p=—j 



(71) 
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This expression is vanishing unless m is even, which implies that j must be an integer. If these conditions arc fulfilled, 
we get, for m = 2m': 2^f 1 = 2j. As for the term m = of the integration over 9 is trivial and restricts j 
to be integral. Gathering these results: 



Zbf = 



" f A f~ TKT 



(t,v) f jew 



It is finally easy to check that X^eN^J + 1) Xj(g) — £so(3)G?)> as desired. 

Let us finally write Zbf as a sum over spin foams, by integrating the group elements g v t- Each of them appears 
four times, once for each face of the tetrahedron t, in the representations, say, ji , j% , j'3 and j'4 . Such a group averaging 
can be written as the identity on the invariant space carrying these representations, Inv(ji, j'2, j'3, ji), under the form 
of a sum over a complete orthogonal basis. The basis elements, called intertwiners, are specified by the choice of a 2 
by 2 pairing and of a (internal or virtual) representation between both pairs. More precisely, we define: 

y ( jl h l )(~iy- m ( 1 ]z k ) (73) 

t-^i \m\ m,2 ml \—m m-$ m^l x ' 



tJ inim2msni4 



which will be shortly denoted i ma {ja, i), although it does not make the pairing explicitly appear. The chosen pairing 
is here j\ with j 2j and the internal representation i. The result of integrating four matrix elements reads: 

f * 

/ dg TT D^ a J (g) = V d 4 i„ la (j a , i) L ria (j a ,i) (74) 
• y su(2) ^ Y 

Then, when an element g v t is integrated out, it leads to two intertwiners respectively due to the 4-simplex v and 
tetrahedron t indices. Since a tetrahedron is shared by two 4-simplices v and v', the previous formula applies 
twice and the tetrahedron is labelled with two intertwiners of representations it v and i tv > carrying indices from the 
tetrahedron ends of g v t and g v n- Moreover, the group element g vt always appears composed with g~, t , so that both 
intertwiners, provided they have the same pairing, are contracted according to the orthogonality relation: 

s x dl 
{m a } 

As for intertwiners with indices from 4-simplices, the five intertwiners corresponding the five tetrahedra meeting at a 
4-simplex v are contracted through the ten faces of v, and thus according to the structure of a 15j-symbol from SU(2) 
recoupling theory [32| . Thus, the partition function can be written [33[: 

z bf= e e n>/ n< n i5 j(j/^o ( 76 ) 

{j/eB}{i*} / * « 



B. Insertions of areas and 3d angles 

The data appearing in the sum over spin foams are usually interpreted through the canonical quantization [34[, 
[Tlj . In particular, the representations jf labelling dual faces stand for lengths of edges in 3d (this is in fact the 
very definition of the Ponzano-Regge model) and areas of triangles in 4d. Using the loop quantization, the spectrum 
of these operators can be computed at the quantum level and their square scales as jf(jf + 1) [33 ] . These results 
can also be obtained from the covariant quantization, by computing quantum average with the functional integral 
and insertions of states exciting the dual faces. It has been properly done recently in 3d (36|. We do not propose 
this treatment here, the issue of the gauge-fixing being more subtle than in 3d, due to the reducibility of the gauge 
transformations [28[ . But we are rather interested in the translation of insertions of local observables into the language 
of spin foams. 

First consider a given dual face /*, and insert into Zbf, (|67|) . the local observable A 2 ^: 



zbf = J n d ^t n n [^/^k^ »/* e-^*)] n d^A^f^ [^(v) a 



{t,v) (f,t) (/,«) jW 



I 2 e *A s 



(77) 
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The result (f?2")) is naturally not modified for / ^ /*. As previously, the variables nf* t are eliminated using the gluing 
relations, except one which leads to (|68|) with g — gf* and — Of*. The representation jf* is also constrained to be 
integral because of the integral over 9. The change comes from the integral over the area: 

dA A 2 J m (A) = m 2 - 1, for m e 2N (78) 

The important point is that, as before, the sum over the angular momentum to/*, which is the variable dual to Of*, 
can be explicitly performed. In addition to the standard dual face weight, 2j + 1, it leads to the insertion into ([72")) 
of the following quantum area: 



.4 



4 

J ~ 3 



3 



i(i + i)-4 =<i(4- 4 ) (79) 



3 



This spectrum gives a negative value for j = 0, would vanish for j — 1/2, and then increases like + !)■ Up to 
a coefficient 1/3, this corresponds to the result obtained in [36j in 3d, with a precise formulation including a gauge 
fixing process. As explained by the authors, the result crucially depends on measure factors and the way curvature is 
described in the lattice theory. It is thus an interesting result, given the different setting chosen here, which involves 
different measures and computations. 

Let us come back to the expansion (|69|). For a given area, all U(l) modes are involved. But after integration over 
A, we are free to interpret each component e m6 ^ 2 as a Regge action with deficit angle and quantized area m. Areas 
are usually rather quantized as spins of SU(2) representations. And indeed, in the previous computations, we see that 
the sum over m can be explicitly performed, leading to standard results. However, in more intricate situations, this 
sum may not be possible to perform and m may be interpreted as a quantized area. To convince oneself, let us apply 
to the same expansion to the action i™ (|55|) : 



ezp(iA f sm(?±6} + jt-9jj) = £ J mf (A f ) e h + ™ f e+ e h-™f<>} ( 80 ) 



Defining = 7±rri/, we obviously have: 



m f = m f h- = m+/ 7+ (81) 



This is precisely the relation usually written between the spins quantizing the areas, [16j, [151 ]. We would have indeed 
written this relation if we have used the bivectors By' = b+J t © b_l t , with bM t = ±7±b±/t to define the theory. 
Diagonal simplicity would have then implied: A_//7_ = A + f/j + for the newly defined self-dual and anti-self-dual 
areas. This strongly suggests that the information concerning the areas are recorded in the spin foam into the U(l) 
modes of the Regge action. 

Let us now study the insertion of the 3d angle, cos which is the angle between the directions of bf t and bf t 

within t. Consider two triangles /i and f 2 in t. This angle can be expressed as: Dqo (n 1 " 1 n 2 ). We naturally expect 
its translation into spin foams to involve the intertwining representation i t , since j\ and j% are coupled through i tl 
if the pairing of the intertwiner is well chosen. We thus write: D$ (n^ 1 n 2 ) = J2a=-i ^oa^i" 1 ) ^Ao( n2 )- After 
integration over areas and dihedral angles, the weight for f\ is: 

{ E d h D^\n^g xni )+ £ £ d h (D^{n^g ini ) + (-l) mi X»^_^(nT V^)) } D^fa 1 ) (82) 



As before, though it has not been explicitly precised, since D_ m _ m {g) = Dm' m (g) — Dm' m (g^ 1 ), the partition function 
sums over the orientation of each dual face with a coefficient (— l) m (independently of the insertion of the 3d angle 
observable, see expression ((70)) ). The result of (p8"|) is now modified and the relevant quantity becomes, for f^: 

L i2) dn ^VV) = - ( -iy- (i \ J) Ec-i^c-i) 1 - 4 (i a I - A ) D<£i9) (83) 

The quantity for / 2 can be obtained using D ( jl(n) = {-l^D^^n- 1 ). Because of the sum over the orientation in 
(fB"2")) , we have to evaluate the properties of (|83"|) under the change / — ► —I. Given the properties of the 3jm-symbols, 
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FIG. 2: The left figure shows the graph weighting the tetrahedron t and resulting from the insertion of the 3d angle between 
/i and /2 in t. In contrast with pure BF theory, it does not simply consist of the contraction of the intertwiners coming from 
the integrations of gi and g2 because of the insertion of a link in the spin 1 representation between the links of /i and fi. This 
graph reduces to the evaluation of a 6j-symbol if the pairing is conveniently chosen, displaying the role of the intertwiners in 
the quantization of the 3d angle. 



one sees that only the sign is changed. In contrast with Zbf, it implies that the variable rrif, dual to the deficit angle 
6 f , has to be odd for fi and , and then that the representations ji and ji are in N + ^ . Notice that the first term 

in 1|82|). corresponding to m — 0, does not contribute since the 3jm-symbols ^ q are zero for j £ N. Like in pure 

BF theory, the sum over the angular momenta mi and m2 can be explicitly performed thanks to the special form of 

the 3jm-symbols. Indeed, £jU 1/2 (-l) J '- fc (ifc { o) = JU+V 

The indices a\, a 2 and b±, 62 of the matrix elements (gi) and D^^j (172) are, after integration over the holonomy 
degrees of freedom, carried by two intertwiners, with different internal representations i and i' taking place at the 
tetrahedron t. In the absence of the observable insertion, one has instead of the 3jm-symbol of (f83|) a Kronecker 
delta, Sab- Now, the contraction of these two intertwiners is not made according to the orthogonality relation. A link 
in the spin 1 representation is inserted between the links j± and J2 of the two faces, as shown in figure O However, 
if the pairing is well chosen, the equality i — i' still holds, and the graph reduces to that of a 6j-symbol. Assuming 
without loss of generality that /1 and f% are similarly oriented and gathering the numerical coefficients, the insertion 
of cost/^j, translates into the insertion into Zbf written as a sum over spin foams of the following quantity: 

(-i) Ji+J2+It / 44 ~ 

16 V + VMh + l) 

where the quantity into brackets denotes the 6j-symbol represented in figure [2] with a tctrahedral graph. The angle 
between two triangles is thus quantized as a 6j-symbol. The latter not only involves j\ and ji which contain the areas 
of fx and /2, but also the representation i t of the intertwiner at which the triangles meet. The intertwining data {it} 
of Z-qf thus encode the correlations between the directions of triangles. 



31 ji 1 
h h it 



3i,h G - 



(84) 



Conclusion 



Precise relations between lattice BF theory, supplemented with the simplicity constraints, and area-angle Regge 
calculus are given. The expressions of the constraints in area-angle Regge calculus and of the (4d) dihedral angles 
in terms of the 3d angles take a simple form through a formulation of the gluing relations and cross-simplicity with 
group elements. The dihedral angles enable to express holonomies in terms of bivectors, and appear as part of a 
u(l) © u(l) discrete connection. We propose two actions (the difference being rather a technical ambiguity) which 
include the Immirzi parameter. It might be interesting to study some other ambiguities, such as choosing a higher 
spin representation instead of the fundamental. 

Our work leads to a mixed setting, combining variables from both Regge calculus and lattice BF gauge theory, and 
whose geometric content is crystal-clear. As emphasized in the introduction, we believe that it is important to have a 
(compactified) Regge-like action for which the introduction of the Immirzi parameter makes sense, in particular while 
aiming at understanding its nature at the quantum level. The Immirzi dependence here disappears from the classical 
equations of motion for non-degenerate configurations. The case of degenerate configurations is less clear, but the 
recent classification results of [191 ] are certainly helpful in this regard. 

Finally, insertions of areas and 3d angles are converted into spin foam language. As expected, areas are quantized 
according to the Casimirs of triangle representations. Their quantization in the presence of the simplicity constraints 
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is briefly discussed. The insertion of a 3d angle corresponds to the insertion of a 6j-symbol, whose data are the spins 
of the two involved triangle representations and the tetrahedral intertwining spin. 

Let us briefly argue that the present setting is well suited to study spin foam models from discretized path integrals. 
We can choose and completely control how cross-simplicity and parallel transport relations are taken into account. 
Indeed, instead of the strong gluing using delta functions over SU(2), {25} , the parallel transport relations can be 
inserted with another weight. It should be peaked around the identity, so that the equations of motion give the correct 
rules, and take the following form: 

j 

Seinjl gwriff e^ 6 ^"') = ^ f jm (e) {j, m\ nj? g tt < n ft . ei e ^ a * \j, m) + c.c. (85) 

j G N m =-j 

Only diagonal elements of the representation matrices are considered, the integrals over the angles 0y^ v and tpft leading 
anyway to this situation [38j |. The strong gluing corresponds to the choice fj m = 2j + 1. e is a generic variable which 
suggests to consider families of models. For instance, we may take S e to be a one-parameter family of Gaussians of 
width e: 

Uk) = ^e^ k3 \ with f jm (e) = W^) ( 86 ) 

where N and N' are normalisations and Ij is a modified Bessel function, defined by Ij(z) — PJj(-iz). Notice that 
this is quite similar to the setting used in semi-classical computations of correlation functions from spin foam models 
[37l ]: one inserts a boundary state which assigns a mean geometry to the boundary of the spacetime region under 
consideration. The above functions precisely assign, as we have shown, mean dihedral angles between the boundary 
tetrahedra for each 4-simplex. The spin foam models resulting from different choices of 5 will all exhibit similar 
structures, since spin foams are based on the integrations of the variables rift and g v t- From this perspective, the 
coefficients fj m may be seen as a measure ambiguity, exactly like that affecting the definition of the EPRL model. 
While the differences between two spin foam models differing by their sets of coefficients fj m may be rather arcane, 
computing the corresponding functions 8 would teach us precisely about these differences. A crucial issue is then 
obviously the ability to resum these coefficients after the integrations. However, some interesting choices of functions 
can prevent us from such a challenge, such as coefficients of the form: fj m oc Sj\ m \, as will be shown in [3a |. In 
any case, the possibility of easily tuning the fluctuations around the gluing relations is a promising tool towards the 
understanding of the content of spin foam models. The same approach can be applied to cross-simplicity, as written 
in the third equation of (|3"U)) . The resulting spin foam models will be presented elsewhere [H[ and in particular the 
way the new spin foam models for quantum gravity arises. 
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